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We calculate the probability to find exactly n eigenvalues in a spectral interval of a large random 
N x N matrix when this interval contains s N eigenvalues on average. The calculations exploit 
an analogy to the problem of finding a two-dimensional charge distribution on the interface of a 
semiconductor heterostructure under the influence of a split gate. 
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Random-matrix theory found its applications in nu- 
merous branches of physics. Among them there are the 
statistical theory of slow neutron resonances, the theory 
of chaotic systems, the properties of small metallic sam- 
ples (mesoscopic physics) etc. [1]. Random matrices were 
proven to be a good approximation for the Hamiltonians 
of such systems. 

The random-matrix theory studies statistical proper- 
ties of spectra of large matrices whose elements have a 
given random distribution. Much interest attracted the 
properties of a spectral interval containing on average s 
eigenvalues, where s is a small fraction of total number 
of eigenvalues N. The actual number n of eigenvalues in 
this interval fluctuates from one realization of a random 
matrix to another. Among different statistical charac- 
teristics of these fluctuations, two particular ones were 
the objects of intensive study. One of them was the vari- 
ance of n: (8n 2 ) = ((n — s) 2 ). Using a thermodynamical 
argument Dyson [2] showed that 
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where [3 = 1,2, or 4 for the three possible ensembles 
of random Hamiltonians: orthogonal, unitary, and sym- 
plectic, respectively. The constant Bp was calculated by 
Dyson and Mehta [3]. It is different for different ensem- 
bles, but is generally of order unity. Here we quoted 
Dyson's result in its simplest form when possible degen- 
eracies and series of non-interacting levels are neglected. 
The fact that (8n 2 ) <C s is the consequence of a so-called 
level repulsion. 

Another popular quantity was the probability density 
P(d) to find two consecutive levels separated by a dis- 
tance of d average level spacings. For the orthogonal en- 
semble P(d) is well approximated by the famous "Wigner 
surmise" P(d) = (ird/2) exp(-7rrf 2 /4) [4] if d is not too 
large. The level repulsion can be seen in the fact that 
P(d) — ► as d — ► 0. For large d the asymptotic behavior 
of P(d) is different from that of "Wigner surmise" [2,1]: 
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In this paper we calculate analytically the asymp- 
totical behavior of a more general and more informa- 
tive quantity. This quantity designated by Mehta [1] as 
Ep(n,s) is the probability to find exactly n eigenvalues 
in a spectral interval containing s eigenvalues on aver- 
age. It is easy to see that P(s) = d 2 Ep(0, s)/ds 2 and 
(hn 2 ) = X^=o( n — s ) 2 Ep{n, s )- The idea behind our cal- 
culation is to push further the analogy of the problem at 
hand to the thermodynamics of a Coulomb gas. It was 
this analogy that Dyson used to obtain Eq. (2). As an 
important result we would like to present the expression 
for Ep(n, s) in the limit \n — s\ <C s: 
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where 6n = n — s and B weakly depends on 6n. This dis- 
tribution is nearly Gaussian for \6n\ <C Ins, and Eq. (1) 
follows. The Gaussian form agrees with the conjecture 
put forward by Al'tshuler et at [5]. For \6n\ ^> Ins the 
probability of a fluctuation 6n in the eigenvalue number 
is, however, significantly smaller than as it would be de- 
termined from a Gaussian distribution with the variance 
given by Eq. (1). 

Let us now turn to the derivation of the general ex- 
pression for Ep{n,s). Suppose we consider the Gaussian 
orthogonal ensemble E\q defined as the ensemble of real 
symmetrical matrices H such that their elements Hij are 
statistically independent and the probability Q(H)dH 
that a system of E\q will belong to the volume element 
dH = Yli<j dHij is invariant under real orthogonal trans- 
formations 
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for any real orthogonal matrix S. The joint probabil- 
ity density function for the eigenvalues of H has the 
form [1] 



Q(xi 



) = C N JJ \xi - xj\exp 



N 



^(axi - b) 2 



(5) 



1 



The quantities a,b, and C'n are some constants. Eq. (5) 
can be rewritten as 
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where 
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and T = 1. Hence, the probability Q is nothing else as 
corresponding Gibbs' factor for a gas of N point charges 
in thermodynamic equilibrium at temperature T. Each 
pair of charges i and j exhibits electrostatic repulsion 
logarithmically dependent on the distance \x{ — Xj\ be- 
tween them. The N charges are confined by an external 
potential <f> c {x) = {ax — b) 2 . Once this analogy is estab- 
lished one can apply the methods of the classical ther- 
modynamics to calculate various statistical quantities of 
random matrix spectra. 

The generalization of Eq. (6) to the other two types of 
ensembles requires only change in the value of the tem- 
perature: T = 1//3. Introduce the partition function 
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The factor 1/N\ corresponds to our treatment of charges 
as indistinguishable particles [2]. The quantity Ep{n,s) 
is the ratio of two values of i\>p calculated for two dif- 
ferent domains of the integration. The first domain is 
determined by the condition that exactly n charges be- 
long to the given interval of length s. The second value 
corresponds to the unrestricted domain of integration. 

When both n and s are large one attempts to aban- 
don the discreet formulation of the problem in favor of 
the continuum one. Following Dyson we make three as- 
sumptions: (i) there is a macroscopic eigenvalue density 
function p{x); (ii) for a given p(x) the free energy of the 
gas is composed of two parts 



F = V 1 + V 2 
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where V\ is the potential energy in the mean-field ap- 
proximation 
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and V'z is the contribution depending only on the local 
density 



V 2 = p(x)fp[p(x)]dx, 
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f[p] being the free energy per particle of a Coulomb 
gas with average density p in equilibrium; (iii) the over- 
whelmingly dominant contribution to the integral (8) 



comes from configurations not deviating significantly 
from an "optimal" density fluctuation, namely that func- 
tion p(x) which makes F a minimum subject to 
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where t = s/2. 

Another simplification may be obtained by noting that 
if s <<^ N then we study only a small fraction of the whole 
spectral interval. In this case the explicit form of the con- 
finement potential in the expression (7) is irrelevant and 
can be chosen according to our needs. We find it the most 
convenient to think that the confinement is due to a com- 
pensating background of unit density. The appropriate 
modification of Eq. (10) is then 
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As Dyson showed [2] the free energy density f[p] is 

f[p] = (l/(]-l/2)p\np. (14) 

The first term is just the entropy multiplied by the tem- 
perature and the second term is the correlation energy 
per unit length. 

The probability Ep{n,s) is expressed as 
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where F is given by Eqs. (9, 11, 13, 14) and the minimum 
is sought in the class of continuous non-negative functions 
p(x) under the condition (12). This optimal fluctuation 
is very close to the one, which ensures the lowest value 
of the mean-field part V\ of the free energy F. However, 
the optimization have to be done taking into account the 
local part V 2 as well, and we will make an error of order 
0(ln|(5n|) as long as s > 1 and \6n\ ^> 1. As discussed 
by Dyson [2] we lost such order terms anyway upon the 
transition from the discreet to the continuum formulation 
of the problem. 

Using the variational principle we find that the elec- 
trostatic potential 
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created by the optimal density fluctuation must satisfy 
the following conditions 
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\x\ < t, p(x) > 
\x\ < t, p(x) = 
\x\ > t 
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if n < s and similarly, 
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(x) = V, \x\ < t, p(x) > 
< V, \x\ > t, p(x) = 
= 0, \x\ > t 
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otherwise. The meaning of these conditions is that the 
electron gas breaks into "metallic" regions, where the 
potential is perfectly screened, and "insulating" regions, 
where the charge density p(x) vanishes and there is no 
screening. 

Consider the case n < s first. We notice that this prob- 
lem is equivalent to another one, which at first glance ap- 
pears completely different. We are referring to a problem 
of finding the charge distribution of a laterally confined 
two-dimensional electron gas (2DEG) in a heterostruc- 
ture. The reason for this is as follows. Due to a large 
aspect ratio in some of such devices, they may be consid- 
ered as translationally invariant in one of the dimensions. 
Then the usual Coulomb interaction of two-dimensional 
charges leads to the logarithmic interaction in terms of 
the one-dimensional charge density. This is precisely the 
type of interaction that we are having in our Coulomb 
gas model. 

Consider a simplified model of the split-gate device 
on the GaAs/Afj.Gai-^As heterostructure studied in 
Ref. [6]. In this model (see Fig. 1), ionized donors and 
the two-dimensional electron gas are characterized by 
continuous charged densities. The donors constitute a 
uniform positive background compensated by the 2DEG. 
The split gate is represented by two semi-infinite metal 
planes separated by the gap of width 2t centered at x = 0. 
The system is translationally invariant in ^-direction and 
all the charges and the gate are in the plane z = 0. 




FIG. 1. The model of the split-gate device. The two shaded areas 
represent the gate in the shape of two semi-infinite planes. Be- 
neath the gate the two-dimensional electron gas located in Z = 
plane is depleted. The density of the electron gas does not vanish 
only in a narrow strip (grey area) between the two halves of the 
gate. 

Suppose a voltage difference V is applied between the 
gate and the 2DEG. The 2DEG will be depleted under 
the gate and near its edges. Therefore, near the edges 
of the gate the charge density is due to the background 
only. If we subtract the total charge density from the 
background charge density, then this new quantity will 
be non-zero in three disconnected regions: in a central 
strip of 2DEG and in the metallic gate. Having in mind 



the analogy discussed above, we designate this quantity 
by the same symbol p(x) as the eigenvalue density of a 
random matrix. Let the width of the strip of the 2DEG 
be 26, then the p(x) in units of the background charge 
density is [6] 



p(x) = Re^/(b 2 -x 2 )/(t 2 -x 2 ). 
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The consequence for the distribution of the eigenvalue 
density in the optimal fluctuation is now that it is 
clearly given by the same formula (19). The param- 
eter 6 < t is to be found from Eq. (12). The graph 
of p(x) is shown in Fig. 2(a). Note that for n = s 
6 = and p(x) = in the interval —t<x<t. 
This particular case was studied by Dyson [2] to 
obtain the asymptotic form of Ep(0,s) and P(s). 
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FIG. 2. The optimal eigenvalue density fluctuation for (a) 71 < S 
(b) Tl > S. 

To calculate the free energy corresponding to the fluc- 
tuation we need to know the distribution of the elec- 
tric field in the Coulomb gas model. It can also be 
obtained knowing the corresponding distribution for the 
split-gate device. Namely, the field in the Coulomb gas 
is by the factor of two weaker. The reason for this is 
that the potential due to real two-dimensional charges 
is —2 f p(y)ln \x — y\, whereas the potential (16) in the 
Coulomb gas model is twice as small. With this correc- 
tion the electric field in the Coulomb gas model is given 
by the expression 



-d<f>(x)/dx = 7T InV(& 2 - x 2 )/(t 2 - x 2 ), 
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FIG. 1. Fogler and Shklovskii "The Probability..." 
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The electric field vanishes everywhere where p(x) is non- 
zero in agreement with Eq. (17). 

The case n > s can be reduced to a similar electrostatic 
problem with the gate in the shape of an infinite strip of 
width 2t to which a positive voltage V is applied. The 
eigenvalue density and the electric field in the Coulomb 
gas are given by the same expressions (19,20) but now 
b > t and the regions of zero p(x) are formed outside the 
interval \x\ < t [see Fig. 2(b)]. 

The calculation of the parameters of the optimal fluc- 
tuations is now straightforward and the answer can be 
given in terms of the complete elliptic integrals [7]: 
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when n < s and 
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when n > s. Here, k 2 = 1 — k' 2 and k' is defined as 
k' = b/t for n < s and k' = t/b for n > s. Also V2 = V 
in both cases. Eq. (15) now reads 

In Ep(n,s) = -iiV x - (1 - /3/2)V + 0(ln \n - s|). (23) 

To use this formula one have to first find the k corre- 
sponding to given n and s and then substitute this value 
into the expressions for V and V\. For example, for n = 
one finds that k = 1 and then V = irt and V\ = 7r 2 t 2 /4. 
One observes that our formulae indeed reproduce the re- 
sult Eq. (2) upon the substitution d = 2t. 

Consider now the case of small fluctuations: \8n\ <C s, 
which corresponds to <C 1- In this limit the expansion 
of Eqs. (21,22) in asymptotic series yields 



Vi 



V 



2\2 



8 



W) 



2 
tk 2 



tk 2 [l 



111 k + 
0(k 2 )] 



0(k 2 ) 



, 4 1 

ln I + 2 



0(k 2 ) 



(24) 



Using Eq. (23) one recovers Eq. (3). As a possible method 
of the numerical check of the latter, we suggest calculat- 
ing the ratio of two subsequent even moments M'2 m +2 and 
M m of E fj{n, s) defined as M m = ~Y^ =0 {n — s) m Epin, s). 
This ratio is expected to be 
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where a = 2/ir 2 /3 and 1 ^ am <C s 2 . In contrast to the 
same kind of ratio for the Gaussian distribution, where 



it is just equal to unity, in our case it decreases with m. 
For instance, at m ~ s/2a it should be about 0.5. 

For completeness we provide also the asymptotic form 
of In Ep(n, s) for large positive 8n, i.e., for n ^> s: 
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One observes that the main terms in both Eq. (3) and 
Eq. (26) are quadratic in 8n and contain large logarith- 
mic factors. This can be explained in terms of the elec- 
trostatic analogy as follows. The quantity In Ep(n,s) 
is up to a constant the work required to charge a two- 
dimensional capacitor by 8n charge units. This work is 
equal to 8n 2 /2C, where C is the capacitance. In the case 
of \bn\ <<^ s the width s of the central plate of this ca- 
pacitor is much larger than the gap between the plates 
and C is logarithmically large. In the case of n ^> s the 
situation with geometrical parameters of the capacitor 
is reversed, and its capacitance is thus inversely propor- 
tional to a large logarithmic factor. 

In conclusion, we studied the probability Ep{n,s) to 
find a given number n of eigenvalues in an interval of a 
random matrix spectrum defined by the condition that 
this interval contains s eigenvalues on average. We cal- 
culated the asymptotical behavior of Ep{n,s) for large 
s. It is found to be Gaussian for small fluctuations of n 
around its average value and to decay faster than Gaus- 
sian at \6n\ > Ins. We suggested a method of the nu- 
merical verification of our result based on calculating of 
large order moments of Ep{n, s). 
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